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The ever-increasing knowledge to the structure of various real- world networks has uncovered their 
complex multi-mechanism-governed evolution processes. Therefore, a better understanding to the 
structure and evolution of these networked complex systems requires us to describe such processes 
in more detailed and realistic manner. In this paper, we introduce a new type of network growth 
rule which comprises of adding and deleting of nodes, and propose an evolving network model to 
investigate the effect of node deleting on network structure. It is found that, with the introduction of 
I node deleting, network structure is significantly transformed. In particular, degree distribution of the 

network undergoes a transition from scale-free to exponential forms as the intensity of node deleting 
increases. At the same time, nontrivial disassortative degree correlation develops spontaneously as a 
£N) ' natural result of network evolution in the model. We also demonstrate that node deleting introduced 

» | , in the model does not destroy the connectedness of a growing network so long as the increasing rate 

O |i of edges is not excessively small. In addition, it is found that node deleting will weaken but not 

^* eliminate the small-world effect of a growing network, and generally it will decrease the clustering 

coefficient in a network. 
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Network structure is of great importance in the topological characterization of complex systems in reality. Actually, 
these networked complex systems have been found to share some common structural characteristics, such as the small- 
world properties, the power-law degree distribution, the degree correlation, and so on [3, 0, In the theoretical 
description of these findings, the Watts-Strogatz (WS) model Q provides a simple way to generate networks with 
the small- world properties. Barabasi and Albert (BA) with a somewhat different aim, proposed an evolving 
network model to explain the origin of power-law degree distribution. In this model, by considering two fundamental 
fS |' mechanisms: growth and preferential attachment (PA), power-law degree distribution emerges naturally from network 
evolution. Based on the framework of BA model, many other mechanisms were introduced into network evolution 
to reproduce some more complex observed network structures [1, 0, H, OH EL EH EH 0, EE El> E3> sucn as the 
degree distribution of broad scale and single scale p , as well as the degree correlation [17| ■ These further studies 
show that real networked systems may undergo a very complex evolution process governed by multiple mechanisms 
on which the occurrence of network structures depends. Therefore, to get a better understanding of the structure and 
evolution of complex networks, describing such processes in more detailed and realistic manner is necessary. 

In the BA's framework, the growing nature of real-world networks is captured by a BA-type growth rule. According 
■ to this rule, one node is added into the network at each time step, intending to mimic the growing process of real 
systems. This rule gives an explicit description to the real-network' growing process which, however, can in fact be 
much more complex. One fact is that in many real growing networks, there are constant adding of new elements, but 
accompanied by permanent removal of old elements (deletion of nodes) [HI, Ell [2(| HH HH, HH] • Take the food webs for 
a example: there are both additions and losses of nodes (species) at ecological and evolutionary time scales by means 
of immigration, emigration, speciation, and extinction [l8|. Likewise, for Internet and the World Wide Web (WWW), 
node-deleting is reported experimentally in spit of their rapid expansion of size [13. 120. |21|. l22t l23j| . In the Internet's 
Autonomous Systems (ASs) map case, a node is an AS and a link is a relationship between two ASs. An AS adding 
means a new Internet Service Provider (ISP) or a large institution with multiple stub networks joins the Internet. An 
AS deleting happens due to the permanent shutdown of the corresponding AS as it is, for example, out of business. 
Investigations of the evolution of real Internet maps from 1997 to 2000 verified such network mechanism [l9|, I20I [2l|. 
The same is for the evolution of WWW, in which the deletions of invalid web pages are also frequently discovered 
[H, HII- In most cases, the deletion of a node is also accompanied with the removal of all edges once attached to it. 
These facts justify the investigation of node-deletion's influence on network structure. In this paper, we introduce a 
new type of network growth rule which comprises of adding and deleting of nodes, and propose an evolving network 
model to investigate the effect of node deleting on the network structure. Before now several authors have proposed 
some models on node removal in networks, such as AJB networks in which a portion nodes are simultaneously removed 
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from the network [241 ) . and also the decaying [25j and mortal [26[ networks, which concerns networks' scaling property 
and critical behavior respectively. Sarshar et al (27j investigated the ad hoc network with node removal, focusing 
on the compensatory process to preserve true scale-free state. They are different from present work, in which node 
deleting is treated as an ubiquitous mechanism accompanied with the evolution of real-world networks. 

This paper is organized as follows. In Section HH an evolving network model taking account of the effect of node 
deleting is introduced which reduces to a generalized BA model when the effect of node deleting vanishes. Then 
the effect of node deleting on network structure are investigated in five aspects: degree distribution (Section 
degree correlation (Section ITV]) . size of giant component (Section [Vj , average distance between nodes (Section [VlJ 
and clustering (Section IVIip . Finally, Section I Villi presents a brief summary. 



II. THE MODEL 



We consider the following model. In the initial state, the network has too isolated nodes. At each time step, either 
a new node is added into the network with probability P a or a randomly chosen old node is deleted from the network 
with probability Pd = X — P a , where P a is an adjustable parameter. When a new node is added to the network, it 
connects to to (to ^ mo) existing node in the network according to the preferential probability introduced in the BA 
model which reads 

n - ka + 1 m 

where k a is the degree of node a. When an old node is deleted from the network, edges once attached to it are 
removed as well. In the model, P a is varied in the range of 0.5 < P a < 1, since in the case of P a ^ 0.5 the network 
can not grow. In order to give a chance for isolated nodes to receive a new edge, we choose preferential probability 
IL a proportional to k a + 1 Q- Note that when P a = 1, our model reduces to a generalized BA model [28| . 

To get a general knowledge to the effect of node deleting on network structure, firstly, a simple analysis to the 
surviving probability D(i,t) is helpful. Here, D(i,t) is defined as the probability that a node is added into the 
network at time step i, and this node (the ith node) has not been deleted until time step t, where t ^ i. Supposing 
that a node-adding event happens at time step i , and the probability that the i'th node has not been deleted until 
time step t is denoted as D'(i' , t). Then, due to the independence of events happened at each time step, it is easy to 
verify that D'(i',t + 1) = D'(i',t)[l - (1 - P a )/N(t)} with D'(i',i') = 1, where N(t) = (2P a - l)t is the number of 
nodes in the network at moment t (in the limit of large t). In the continuous limit, we obtain 

dD'{l',t) (I -P a ) 

—at~ = -J2p—^)i 1 (,) ' () 

which yields 

£ N -(l-P a )/(2P a -l) 



D'(i',t)=[-\ . (3) 



Thus to get the D(i,t) we should multiply D'(i',t) with P a , i.e. 

/ f \ -(1-Pa)/(2P„-1) 

£>(*,*)= Pal 7 ■ (4) 



One can easily find that D(i, t) decreases rapidly as t increases and/or as i decreases provided 0.5 < P a < 1. It is well 
known that highly connected nodes, or hubs, play very important roles in the structural and functional properties of 
growing networks [l], 0, H|- The formation of hubs needs a long time to gain a large number of connections. As a 
consequence, according to Eq. ([5]), a large portion of potential hubs are deleted during the network evolution. Thus 
it can be expected that the introduction of node deleting has nontrivial effects on network structure. In the following 
we show how network structure can be effected by the node deleting introduced in present model. 



III. DEGREE DISTRIBUTION 



The degree distribution p(k), which gives the probability that a node in the network possesses k edges, is a very 
important quantity to characterize network structure. In fact, p{k) has been suggested to be used as the first criteria 
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to classify real- world networks Q. Therefore it is necessary to investigate the effect of node deleting on the degree 
distribution of networks firstly. Now we adopt the continuous approach [2^ | to give a qualitative analysis of p(k) 
for our model with slight node deletion (i.e., when Pd is very small). Supposing that there is a node added into the 
network at time step i', and this node is still in the network at time t, let k(i',t) be the degree of the z'th node at 
time t, where t ^ i' . Then the increasing rate of k(i',t) is 



dk(i',t) 

at 



P„m 



k(i',t) + 1 
S(t) 



(1 



where 



M',t) 
N(t) • 



S(t)=J2D'(i',t)[k(i',t) + l] 



(5) 



(6) 



and the denotes the sum of all i' during the time step between and t. It is easy to verify that the first term 
in Eq. ([5]) is the increasing number of links of the i'th node due to the preferential attachment made by the newly 
added node. The second term in Eq. © accounts for the losing of a link of the i'th node during the process of node 
deletion, which happened with the probability k(i',i)/N(t). 
Firstly we solve for the S(t) and get 



S(t) = (2P a -l)(2P a m+l)t 
(see the Appendix for details). Inserting Eq. |(7J) back into Eq. ([5]), one gets 

dk(i',t) _ Ak(i',t) + B 



dt 



where 



and 



A = 



2P a 2 m ~P a m + P a -l 
(2P„-l)(2P«m + l) 



B = 



Pam 

(2P Q -l)(2P a m+l)' 



When Ak + B > 0, the solution of Eq. © is 

k(i',t) = 



(Am + B) - - B 



(7) 



(8) 



(9) 



(10) 



(11) 



Now, to get the probability p(k, t) that a randomly selected node at time t will have degree k, we need to calculate the 
expected number of nodes Nk(t) with degree k at time t. Then the p(k, t) can be obtained from p(k, t) = Nk(t)/N(t), 
where N(t) is the total number of nodes at time t. Let Ik{t) represent the set of all possible nodes with degree k at 
time t, then one gets 



p(k,t) 



1 



Nk(t) 

N(t) N(t) 



(12) 



iei k (t) 



In the continuous-time approach, the number of nodes in Ik(t) is the number of i's for which k ^ k(i,t) ^ k + 1, 
and it is approximated to \dk(i,t)/di\^_ i , where it is the solution of the equation k(i,t) = k. To proceed with our 
analysis, now we make the approximation that all nodes in Ik(t) have the same surviving probability D(if.,t) [44| . 
Under this mean-field approximation, Eq. (|12j) can be written as 



p(k,t) 



N(t) 



D(i k ,t) 



dk(i,t) 



From Eq. (|11| . we obtain 



Ak + B 
Am + B 



di 



-l/A 



(13) 



(14) 
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FIG. 1: P™ in [defined in Eq. ((20j)] as a function of m. 



then 

-1 



dk(i,t) 



di 

Inserting Eq. (jT^J) back into Eq. (g]) we get 



(Am + B) 1/A t(Ak + By {A+1)/A . (15) 



Inserting Eqs. ([T5j) and HU) into Eq. (H]), and noting that N(t) = (2P Q - l)t, we get 

p{ki t ) = -JiL_ (Am + B ) {B - A+1)/A (Ak + B )- {B+1)/A , (17) 

which is a generalized power-law form with the exponent 

7= ^±1 = 2+ Pam + 1 . (18) 

' A 2P a 2 m - P a m + P a -l 

We point out again that equation (fTTj) is only valid when Ak + B > 0, which translates into A > 0, i.e. 

2P Q 2 TO _ P a m + P a - 1 > 0. (19) 

Considering that P a > 0.5, Eq. (fT9j) is satisfied when 



p . pmin _ (m - 1) + \/7n 2 + 6m + 1 

In Fig.[Tl we plot p™ n as a function of m. One can see from Fig.Q]that the curve divides our model into two regimes. 
(i) P a > P™ m : in this case Ak + B > and equation (fTTj) is valid. Thus, the degree distribution of the network 
p(k) exhibits a generalized power-law form, (ii) P a > p™ ln \ In this case Ak + B > can not be always satisfied 
and equation (jlip is not valid. Therefore, our continuous approach fails to predict the behavior of p(k), and we will 
investigate it with numerical simulations. The P™ m (m), as one can find from Fig. [TJ decreases with the increase of 

771. 

In the power-law regime [P a > P™ m (m)], the behavior of p(k) is predicted by Eqs. (fT7|) and (fT5|) , which are obtained 
using a mean- field approximation [Eq. (|13[) ] . One can easily verify that such approximation is only exact when P a = 1 , 
in which case Eq. (|18p turns into 7 = 3+ I/ttt,, in good agreement with the results obtained from generalized BA 
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FIG. 2: Cumulative degree distribution P(k) for networks with system size N — 100000 and different values of P a , in logarithmic 
scales. The dash line is power-law fit for P a = 1. The solid line is the exponential fit for P a = 0.51. In the simulation, we 
set mo — m — 5 and each distribution is based on 10 independent realizations. Inset plots the power-law exponential 7 as a 
function of P a . The continuous curve is according to the analytic result of Eq. (|18jl . and circles to the simulation results. 



model studied in Ref [28| . If P a nm (m) < P a < 1, Eqs. (fl7|) and (fTSf still give qualitative predictions for the model: 
with slight node deletion, p(k) of the network is still power-law, and the exponential 7 increases with the decrease of 
P a (inset of Fig. El). 

In remaining regime [P a < P™ m (m)], the limiting case is P a — » 0.5, in which the growth of network is suppressed (a 
very slowly growing one). Similar non-growing networks have been studied, for example, for the Model B in Ref 30], 
and the degree distribution has the exponential form. Here we conjecture that, in this regime, p{k) of our model 
crossovers to an exponential form, which is verified by the numerical simulation results below. 

Now we verify the above analysis with numerical simulations. In Fig. El we give the cumulative degree distributions 
P{k) [3J] of the networks with different P a . As P a gradually decreases from 1 to 0.5, Fig. El shows an interesting 
transition process which can be roughly divided into three stages. (1) 0.9 ^ P a ^ 1: In this stage, the model works in 
the power-law regime and the power-law exponent 7 increases as P a decreases. Inset of Fig. El gives the comparison 
between the value of 7 predicted by Eq. (fT8|) and the one obtained from numerical simulations. One sees that the 
theory and the simulation results are in perfect agreement for P a = 1. As P a decreases, however, the agreement is only 
qualitative and the deviation between theory and simulation becomes more and more obvious. As we have mentioned 
above, such increasing deviation is due to the mean-field approximation used in the analysis. These results tell us that 
slight node deletion does not cause deviation of the network from scale-free state, but only increases its power-law 
exponent. Such robustness of power-low p(k) revealed here gives an explanation to the ubiquity of scale-free networks 
in reality. It should be noted that a very similar robustness has also been found in the study of network resilience, 
where simultaneously deleting of a portion of nodes was taken into account in static scale-free networks [Hj]. (2) 
0.5 < P a ^ 0.6: In this stage, the model works in the regime of P a < P a Tlln (m). As one sees from Fig. El P(k) of 
the network behaviors exponentially. This result indicates that with manifest node deletion, the network will deviate 
from scale-free state and become exponential. (3) 0.6 < P a < 0.9: In this stage, a crossover of the model from the 
power-law regime to the exponential regime is found, in which the P{k) is no longer pure scale-free but truncated by 
an exponential tail. As one can see, the truncation in P(k) increases as P a decreases. 

Besides the power-law degree distribution, it is now known that p(k) in real world may deviate from a pure power- 
law form [l8|, [31], [H, Hi| . According to the extent of deviation, p(k) of real systems has been classified into three 
groups @|: scale- free (pure power-law), broad scale (power-law with a truncation), and single scale (exponential). 
Many mechanisms, such as aging [1, d, H[ , cost @, and information filtering jlQ] . have been introduced into network 
growth to explain these distributions. Here, the results of Fig. El indicate that a modified version of growth rule can 
lead to all the three kinds of p(k) in reality, and it provides another explanation for the origin of the diversity of 
degree distribution in real-world: such diversity may be a natural result of network growth. 
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FIG. 3: Assortativity coefficient r plotted with network size AT, for different P a in the model. In the simulation, mo = m = 5. 
Result of each curve is based on 10 independent realizations. 



IV. DEGREE CORRELATION 



It has been recently realized t hat, besides the degree distribution, structure of real networks are also characterized 
by degree correlations [H, HE H3, HI] • This translates into the fact that degrees at the end of any given edge in real 
networks are not usually independent, but are correlated with one another, either positively or negatively. A network 
in which the degrees of adjacent nodes are positively (negatively) correlated is said to show assortative (disassortative) 
mixing by degree. An interesting observation emerging from the comparing of real networks of different types is that 
most social networks appear to be assortatively mixed, whereas most technological and biological networks appear to 
be disassortative. The level of degree correlation can be quantified by the assortativity coefficient r lying in the range 
— 1 ^ r ^ 1, which can be written as 

M- 1 E i ^-[M- 1 E4fe + fc 1 )] 2 (21) 

for practical evaluation on an observed network, where ji, ki are the degrees of the vertices at the ends of the ith edge, 
with i = 1, . . . ,M [35]. This formula gives r > 0(r < 0) when the corresponding network is positively (negatively) 
correlated, and r — when there is no correlation [45j . 

Recently, Maslov et al [1^] and Park et al [3(3] have proposed a possible explanation for the origin of such correlation. 
They show for a network the restriction that there is at most one edge between any pair of nodes induces negative 
degree correlations. This restriction seems to be an universal mechanism (indeed, there is no double edges in most 
real networks), therefore, the authors of Ref. (4(| conjecture that disassortativity by degree is the normal state of 
affairs for a network. Although only a part of the measured correlation can be explained in the way of Ref. [4(| , this 
universal mechanism does give a promising explanation for the origin of degree correlation observed in real networks 
of various types. 

It will be of great interest to discuss the effect of node deleting on degree correlation. In Fig. [3j we give the 
assortativity coefficient r as a function of network size N, for different P a in our model, for m = 5. As one sees 
from Fig. [3J for each value of P a , after a transitory period with finite-size effect, each r of networks tends to reach 
a steady value. When P a = 1, r — > as N becomes large. This result indicates that networks in the BA model are 
uncorrelated, in agreement with results obtained in previous studies [3o1 . l38j . When P a < 1, nontrivial negative degree 
correlations spontaneously develop as networks evolve. One can see from Fig.[3]that the steady value of r in the model 
decreases with the decreasing P a . In particular, when P a ^ 0.6, the value of r is about —0.1. These results indicate 
that node deleting leads to disassortative mixing by degree in evolving networks. To make such relation more clear, 
in Fig. [4} we plot r of networks in our model as a function of P a , for different m. As the Fig. [3J indicates, when the 
network size is larger than 40000, the assortativity coefficient r is nearly stable. So all results in Fig. 0]are obtained 
from networks with N = 40000. Fig. Q] gives us the same relation between r and P a shown in Fig. [3] What is more, 
it tells us that for a given P a , r will increase with the increasing m. The increment gets its maximum between m = 1 
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FIG. 4: Assortativity coefficient r as a function of P a , for different m in the model. In the simulation, N = 40000. Result of 
each curve is based on 10 independent realizations. 
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FIG. 5: Assortativity coefficient r plotted with network size N, for different P a in the randomly growing network model. In 
the simulation, mo = m = 5 and each curve is based on 10 independent realizations. 



and other values. We point out that this is because when m = 1, the network has been broke up into small separate 
components (see the following section). We can also find from Fig. 0] that the gap between different curves decreases 
with the increasing m and the curves tend to merge at large m. 

Now we give some explanations to the above observations. In the BA model, the network being uncorrelated is 
the result of a competition between two factors: the growth and the preferential attachment (PA). On the one hand, 
networks with pure growth is positively correlated. This is because the older nodes, also tending to be higher degree 
ones, have a higher probability of being connected to one another, since they coexisted earlier. In Fig. [51 we compute 
the assortativity coefficient r of a randomly growing network, which grows by the growth rule of BA-type, while the 
newly added nodes connect to randomly chosen existing ones. As one can see from Fig. [5] that pure growth leads to 
positive r. On the other hand, the introduction of PA makes the connection between nodes tend to be negatively 
correlated, since newly added nodes (usually low degree ones) prefer to connect to highly connected ones. Then 
degree correlation characteristic of the BA model is determined by this two factors. In Fig. O we plot the average 
degree of the nearest neighbor < k > nn as a function of k in the BA model. It is found that nodes with large k show 
no obvious biases in their connections. But there is a short disassortative mixing region when k is relatively small 
(also reported in Ref. [4l|, see Fig. la therein). Such phenomenon can be explained by the effect of these two factor: 
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FIG. 6: Average degree of the nearest neighbor as a function of k for the BA model. In the simulation, N = 10000 and 
m = mo = 5. Result of each curve is based on 1000 independent realizations. 
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FIG. 7: The relative size of the largest component S as a function of P a for m — 2,3, 4, 5. Inset gives the same curve for m — 1. 
In the simulations, N = 100000. All results are based on 10 independent realizations. 



Growth together with PA makes nodes with large k equally connect to both large and small degree nodes, and the 
latter makes nodes with small degree be disassortatively connected. Now, we introduce node-deletion. According to 
Eq. (01, depression of the growth of large-degree nodes also decreases the connections between them, therefore makes 
the correlation negative. We also investigate the effect of node deleting on the r of the randomly growing network, 
and obtained similar results. As one sees from Fig.[5J depression of connections between higher degree nodes causes 
the network less positively correlated, and with stronger node-deletion, negatively correlated. Finally, with regard to 
the effect of in in this relation (Fig. [J}, larger m means more edges are established according to the PA probability 
Eq. fl}. We conjecture that the orderliness of newly added nodes connecting to large degree nodes will be weakened 
by the increasing randomness as m becomes larger, thus leading to a less negative correlation. Such randomness can 
not always increase and, as we see from Fig. 2J for large m, e.g., m > 14, the curves tend to merge together. 
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V. SIZE OF GIANT COMPONENT 



In a network, a set of connected nodes forms a component. If the relative size of the largest component S in a 
network approaches a nonzero value when the network is grown to infinite size, this component is called the giant 
component of the network [3, 0, Q ■ In most previously studied growing models [J H El , due to the BA-type growth 
rule they adopted, there is only one huge component in the network, i.e., 5=1. In this extreme case the network 
gains a perfect connectedness. The opposite case of S = 1 is the extreme of S — 0, in which case the network, 
made up of small components, exhibits no connectedness. Experiments indicate that some real networks seem to 
lie in somewhere between these two extreme: they contain a giant component as well as many separate components 
@, H, 51 El- For example, According to Ref.Q, in May of 1999, the entire WWW, containing 203 x 10 6 pages, 
consisted of a giant component of 186 x 10 6 pages and the disconnected components (DC) of about 17 x 10 6 pages. In 
general, the introduction of node deletion in our model will cause the emergence of separate components even isolated 
nodes in the network. What we interest here is the connectedness of the network. In Fig. [7] we plot the relative size 
of the largest component S in the model, as a function of P a , for m = 2,3,4,5, where to is the number of edges 
generated with the adding of a new node. One sees from Fig. [7] that for any 0.5 < P a < 1, a giant component can be 
observed in the model if to > 1. In addition, for the same P a , S increase as the increase of to. While when to = 1, the 
network is found to be broke up into separate components if P a < 1. For example, when P a = 0.9, S of the network 
with N = 100000 rapidly drops to 0.034. Inset of Fig. [JJ gives the S Vs P a curve for m = 1. These results indicate 
that node deleting does not destroy the connectedness of a growing network so long as the increasing rate of edges is 
not excessively small. 



VI. AVERAGE DISTANCE BETWEEN NODES 



Now we study the effect of node deletion on networks' average distance L between nodes. Here the distance between 
any two nodes is defined as the number of edges along the shortest path connecting them. It has been revealed that, 
despite their often large size, most real networks present a relatively short L, showing the so-called small- world effect 
[H, H, 3, 4J. Such an effect has a more precise meaning: networks are said to show the small- world effect if the value 
of L scales logarithmically or slower with network size for fixed mean degree. This logarithmic scaling can be proved 
for a variety of network models [J H [![ . As we have demonstrated in Section [V] node deleting does not destroy the 
connectedness of the network in our model for any m > 1, since there is always a giant component exists. Here in 
our simulation, we calculate L of the giant component of the network in our model using the burning algorithm [3[. 
In Fig. [51 we plot I as a function of network size N, for different P a in our model. As one can see from the figure, 
for any 0.5 < P a < 1, a logarithmic scaling L ~ IniV is obtained, while the proportional coefficient increases with 
the decrease of P a . Furthermore, for a given N, L increases with the decrease of P a . These results tell us that node 
deleting will weaken but not eliminate the small-world effect of a growing network. 



VII. CLUSTERING 



Finally, we investigate the effect of node deletion on network's cluster coefficient C, which is defined as the average 
probability that two nodes connected to a same other node are also connected. For a selected node i with degree ki 
in the network, if there are Ei edges among its ki nearest neighbors, the cluster coefficient C, of node i is defined as 

Ci = h^liy (22) 

Then the clustering coefficient of the whole network is the average of all individual d. In Fig. [HI we plot C of the 
giant component in the network as a function of network size N, for different P a . As one sees from Fig.|9l for each P a , 
the clustering coefficient C of our model decreases with the network size, following approximately a power law form. 
Such size-dependent property of C is shared by many growing network model [l|, |2|, [3j. Moreover, as Fig. [5] shows, 
for the same network-size N, C decreases as P a decreases. The results of Fig. |H] indicate that node deleting weakens 
network's clustering. 



VIII. CONCLUSION 



In summary, we have introduced a new type of network growth rule which comprises of adding and deleting of nodes, 
and proposed an evolving network model to investigate effects of node deleting on network structure. It has been 
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FIG. 8: Average distance L of the giant component in the network as a function of network size N, for different P a in the 
model. The chose of some parameters: mo = m = 5. These curves are results of 10 independent realizations. 




FIG. 9: Cluster coefficient C of the giant component in the network as a function of network size N, for different P a . In the 
simulation we set mo = m = 5. These curves are results of 10 independent realizations. 



found that, with the introduction of node deleting, network structure was significantly transformed. In particular, 
degree distribution of the network undergoes a transition from scale-free to exponential forms as the intensity of node 
deleting increased. At the same time, nontrivial disassortative degree correlation spontaneously develops as a natural 
result of network evolution in the model. We also have demonstrated that node deleting introduced in our model does 
not destroy the connectedness of a growing network so long as the increasing rate of edge is not excessively small. In 
addition, it has been observed that node deleting will weaken but not eliminate the small-world effect of a growing 
network. Finally, we have found that generally node deleting will decrease the clustering coefficient in a network. 
These nontrivial effects justify further studies of the effect of node deleting on network function which include 
topics such as percolation, information and disease transportation, error and attack tolerance, and so on. 
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APPENDIX: THE CALCULATION OF S(T) 
To get S(t), we multiply both sides of Eq. ([5]) by D'(i',t) and sum up all i' between and t: 

i 

To get the above equation we have used the definition of S(t) [Eq. ©] and the following equation: 

Vfl'(i,()= ( D(i,t)di. (A.2) 
Jo 

The left-hand side of Eq. (|A.1[) can be simplified as: 

^ 9{[fcy,t) + i] p'(i',t)} ^ , ^ diy{i',t) 

= | |E *) + !] ^ f ) | - ^ *) + !] *) 

i 

Substituting the above expression in Eq. (jA.ip , and noting that k(t, t) — m and D(t, i) — P a , we get 

dt (2P a - l)t W 

The solution to the above equation is 

S(t) = {2P a -l)(2P a m + l)t. 
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